epl draft 



Threefold onset of vortex loops in superconductors with a mag- 

^netic core 

o : 

^Mauro M. Doria 1 ' 2 Antonio R. de C. Romaguera 1 - 2 , M. V. Milosevic 1,3 and F. M. Peeters 1 

/"■s Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 

Instituto de Fisica, Universidade Federal do Rio de Janeiro, CP. 68528, 21941-972, Rio de Janeiro RJ, Brazil 
V~) 3 Department of Physics, University of Bath, Claverton Down, BA2 7 AY Bath, United Kingdom 



O 

o 
u 

Oh! 



o 
o 
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PACS 74. 25. Op - Mixed state, critical fields, and surface sheath 
PACS 74 . 78 . Na - Mesoscopic and nanoscale systems 

Abstract. - A magnetic inclusion inside a superconductor gives rise to a fascinating complex 
of vortex loops. Our calculations, done in the framework of the Ginzburg-Landau theory, reveal 
that loops always nucleate in triplets around the magnetic core. In a mesoscopic superconducting 
sphere, the final superconducting state is characterized by those confined vortex loops and the 
ones that eventually spring to the surface of the sphere, evolving into vortex pairs piercing through 
the sample surface. 



> 

r-- 
o 

<N 

r— I 

r-- 



Examples of one- or multi-dimensional collections of 
topological defects arise in various physical systems, such 
as vortices in quantum fluids [1], dislocations in solids 
[2], global cosmic strings [3,4], and polymer chains [5]. 
Changes of symmetry in their structural and dynamical 
properties often lead to very rich phase transitions. In su- 
perconductivity, the structure and symmetry of the vortex 
lines has been a major issue even in the early days, when 
Abrikosov [6] proposed a square pattern for the lattice of 
penetrating flux tubes in bulk superconductors; later on, 
Kleiner et at. [7] found that the most efficient way to pack 
interacting parallel rods is actually through a hexagonal 
pattern. However, in systems with reduced dimension- 
ality, the Abrikosov symmetry of the vortex state can be 
strongly affected by the sample boundaries. In mesoscopic 
superconductors, where the volume to surface area ratio is 
small [8], strong lateral confinement may impose the for- 
mation of multiquanta ('giant') vortices [9,10], while the 
shape of the boundary directs the symmetry of the final 
vortex configuration [11]. 

In this Letter, we consider a novel vortex structure that 
stems from the specific inhomogeneous magnetic field of a 
point like magnetic dipole with moment \i in the center of 
a superconducting sphere. This magnet is the only source 
and sinkhole of vortex lines and sets an axis around which 
the vortex state exhibits a discrete symmetry. Because 
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the sphere does not break this degeneracy many of the 
present results also apply for a magnetic inclusion in a 
bulk superconductor. 

Therefore, the present analysis is very different from ear- 
lier studies of vortex and antivortex structures in super- 
conductor/magnet hybrids [12-14] where magnetic dots 
were placed near, but outside a thin superconducting sam- 
ple. While in the latter case vortices were two-dimensional 
(2D) structures characterized by a coin-like core, here they 
are lines in three-dimensional (3D) space, set on a formal- 
ism that comprises the full 3D Ginzburg-Landau theory, 
thus far more complex than previous works. Consequently, 
novel vortex phases with profound 3D features are found 
for these sets of curved vortex lines, made of confined vor- 
tex loops (CVLs) and external vortex pairs (EVPs). 

For small /i, as shown later, we find that the vortex state 
is made of exactly three confined vortex loops around the 
dipole axis. For an increasing magnetic moment fi more 
elaborate vortex arrangements arise because of the inho- 
mogeneity of the field. This threefold growth of vortex 
loops in our system is shown here to be a consequence 
of the energetic balance. Threefold symmetry has been a 
source of inspiration since the sixties when its discovery 
in hadronic physics has led to quarks as the basic con- 
stituents of matter. Recently, experimental evidence was 
found for the onset of trimer bound states in ultracold 
gases of caesium atoms [15] and Z3 symmetry was also 
proposed as a way to solve baryon number violation in 
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grand unified theories [16]. 

Theoretical concept. — The magnetic field of the 
dipole B = [3(/x • f) — /x]/r 3 produces a mixed state whose 
size can be estimated by simple arguments. Take its mag- 
netic field, at the equatorial plane (/x-f = 0), and compare 
it to the known upper and lower critical field expressions 
obtained for a lattice of parallel flux lines, namely, H C 2 = 
<£>o/27r£ 2 , and H c \ = H C 2 ln/^/2^ 2 , respectively, where £ is 
the coherence length and n is the Ginzburg-Landau con- 
stant. The critical field regions are reached at distances 
r(H c2 ) = (/i//i ) 1/3 £ and r(H cl ) = r(H c2 )(2n 2 / Inn) 1 ^, 
respectively. The condition r(H c2 ) > £ holds for a mag- 
netic moment larger than /io = <3>o£/27r. In this case the 
volume around the magnetic moment, where supercon- 
ductivity is destroyed has a radius larger than the mini- 
mal length scale £, the Cooper-pair size. The feasibility 
of such a magnetic domain with moment /io is verified by 
considering a finite volume instead of a point-like mag- 
netic moment. Let us assume a spherically magnetized 
core with radius vm , and compare its magnetization M to 
the saturation magnetization, defined here to be one Bohr 
magneton per atom, Mo = /jLb/(^o,q/3) = 1.49 tesla, 
where ao « 0.05 nm is the atomic radius. The Bohr mag- 
neton is fiB = ^o^c/^Tr, and r c is the electron classical 
radius. For instance, for £ w 3.0 nm, /io has nearly a mil- 
lion oriented Bohr mag netons (f/r c « 10 6 , £/a « 60). 
In this case the magnetized core must have radius slightly 
larger than the coherence length, namely, tm > l-7£ in 
order to have a magnetization smaller than the saturation 
value (M < Mo). Here, we consider a type II superconduc- 
tor of finite size, such that the London penetration length 
is much larger than the sample and therefore, r(H c i) is 
never reached. The magnetic field streamlines around the 
magnetic moment are described in spherical coordinates 
by r = a sin 2 0, a being a parameter associated to the 
loop size, and 6 the angle with respect to the magnetic 
moment axis. This is a special case of the so-called rose 
curve, firstly described by the seventeenth century Italian 
mathematician, Guido Grandi. 

To investigate the superconducting state of a sample 
with volume V, we minimize with respect to the order 
parameter ip the Ginzburg-Landau (GL) free energy 

^ = /f (lHv-^)^| a -H 2 + ^l 4 ), (i) 

expressed in units of the critical field energy density, 
Fo = H 2 /4tt. Note that in Eq. (1) all distances are scaled 
to £, and magnetic moment /i is given in units of /io- The 
order parameter is kept equal to zero in the center of the 
sphere where the magnetic moment is located, a condition 
that does not affect our results since the previously de- 
scribed H C 2 normal core has radius larger than £ around 
the dipole. Three different and independent numerical 
methods were applied to this GL theory: simulated an- 
nealing minimization [17], integration of the GL equa- 
tions [9], and kinetic energy eigenfunction expansion [18]. 



Fig. I: (Colour on line) Isoplots of the Cooper-pair density 
for a magnetic moment /i//io equal to 15, 20, and 25 (left to 
right) in a superconducting sphere of radius 15. 0£. Red and 
grey surfaces are parts of the same isosurface drawn at 20 % 
of the maximum \ip\ 2 . 

Therefore, the present results were checked to be method 
independent. In what follows, we consider a sphere with 
radius R = 15.0£, and determine the vortex state as the 
magnetic moment /i is increased. 

The threefold symmetry. — As the magnetic dipole 
strongly suppresses superconductivity inside the H c2 core, 
vortex loops nucleate as 3D objects in its vicinity, as 
shown in fig. 1. The following animations show a three- 
dimensional view of the loop states upon rotation of the 
sphere for the /i//io equal to 19 and 28 states, which cor- 
respond to embryonic and fully developed loop states, re- 
spectively: M19.mov and M28.mov (see multLmedia file 
for more details about the animations). Distinctively from 
2D vortices the length of the vortex line plays a funda- 
mental role here. We find the remarkable result - a single 
loop or a pair of loops is not a stable state; three vortex 
loops are energetically favorable. With increasing dipole 
strength vortex loops nucleate in triplets from the H c2 core, 
in a continuous manner, i.e. through a second order tran- 
sition. We should emphasize that such a behavior was 
found regardless of the sample size (for R > 10£) and 
symmetry (e.g. sphere, cube), and is therefore a feature 
that holds also for bulk samples. The threefold growth 
of vortex loops is a consequence of the delicate energetic 
balance among the three-dimensional curved vortex lines. 
We used two methods to check its stability. 

The first method is the kinetic energy eigenfunction ex- 
pansion [18], a very appropriate method to prove this be- 
cause near to the core the order parameter must be small. 
We expanded the order parameter, \I/ = J2 k =i Ck^k in the 
orthonormal eigenfunctions of the 3D kinetic energy oper- 
ator (ih\/sD — e*A/c) 2 ^/2m* = e^k- In cylindrically 
symmetric structures, these eigenfunctions have the form 
ipk=(L,n) — ex P(^0)/n( r 5 #)? where L is the angular mo- 
mentum, (j) is the azimuthal angle around the dipole axis, 
and the index n counts different states with the same L. 
The eigenfunctions f n are real and their corresponding 
eigenvalues are obtained numerically for each L, using the 
Housholder technique. Neumann boundary condition are 
assumed on the surface of the sample. The typical number 
of considered eigenfunctions is £7=10-50, with L G (—5, 5) 
and n G (1,5). To search for the free energy minimum 
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Fig. 2: (Colour on line) The free energy as a function of 
the dipole moment, obtained by the eigenfunction expansion 
method. Insets show the phase plots (blue/red - 0/2tt) of the 
corresponding vortex states, taken in the central cross-section 
perpendicular to the dipole. 



we map the superconducting state into a 2D cluster of 
£1 classical particles (xk,yk) — (^(Cfc)? S(Cfe)), which is 
governed by the energy functional obtained when the ex- 
pression for \I> is substituted in Eq. (1). In other words, the 
energy minimization with respect to the complex variables 
Ck leads to the stable vortex states. Note that the total 
angular momentum of the vortex loops equals zero, thus 
both eigenfunctions with positive and negative L must be 
included [L G (— L maxi L max )]. Maximal considered vor- 
ticity L max determines the maximal achievable number 
of vortex loops. For L max = 1 and L max = 2, we were 
able to stabilize one and two loops respectively, but their 
energy was always significantly higher than the Meissner 
state, as shown in fig. 2. The encircling discontinuity of 
2tt (blue to red/red to blue) in the phase contour plots of 
fig. 2 signals a vortex/antivortex. However, when eigen- 
functions with higher angular momenta were included in 
the calculation, the energy minima associated with just 
one and two loops ceased to exist. Instead, stable states 
with one pronounced loop and two 'satellite' ones (and 
vice versa) were found, yet with higher energy than the 
state with perfect threefold symmetry. 

The second method used to support the stability of 
the three loop state is based on the minimization of the 
Ginzburg-Landau free energy through the simulated an- 
nealing procedure [17] starting from an initial constructed 
state, ^at(<2, (/>), which depends on the following variables: 
the loop size, a, the azimuthal angle around the dipole 
axis, and the number of zeros of this wave function, 
N. This initial state is built such that the condition 
^aK^o, </>o) = holds for a single loop position, a = ao, 
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Fig. 3: (Colour on line) The free energy (arbitrary units) of the 
final relaxed state is shown versus the ratio loop size to radius, 
ao/R, for the initial state. The magnetic moment is equal to 
ji/jio — 25 and the initial states are made up of either one, 
two or three symmetrically arranged loops. If the initial state 
is near to the edge, the onset of external vortex pairs becomes 
possible otherwise not. 



and for a special azimuthal configuration determined by 
the condition exp(ziV(/>o) = 1. Hence this initial state van- 
ishes along TV symmetric loops. Obviously the above re- 
quirements don't uniquely define ^jy (a ,0), but we don't 
discuss the specific wavefunction taken here as we found 
invar iance in our results with respect to this choice. For a 
fixed value of fi the numerical search for the minimum is 
carried and the final state free energy is shown in fig. 3 ver- 
sus the initial loop size, ao/R. Three distinct initial states 
were taken, i.e, TV = 1, 2 and 3 CVLs. Indeed the nu- 
merical procedure yielded just one single final state with 
3 CVLs regardless of the number of CVLs in the initial 
state, thus proving its stability. However this holds only 
for a small ao/R ratio, according to fig. 3. For a large 
ao/R ratio several final states are reached. The proxim- 
ity to the boundaries explains this puzzling fact that for 
ao/R > 0.73 the free energy minimum depends on the ini- 
tial state, as shown in fig. 3. The proximity of N loops to 
the surface of the sphere causes the onset of EVPs during 
the annealing procedure. In fact for the three highest val- 
ues of ao / R, the number of EVPs is exactly the number 
of CVLs in the initial state (TV). This striking correla- 
tion leaves no doubt that an initial state with a CVL very 
near to the edge most easily turns into an EVP during 
the minimization procedure. Regardless of the number 
of EVP all the observed final vortex states of fig. 3 dis- 
play three embryonic confined loops very near to the H C 2 
core, as shown in fig. 4. The stability of states with other 
than three loops is only possible by explicitly breaking the 
azimuthal symmetry, such as in case of a rigid lattice of 
magnetic moments inside a bulk superconductor [19,20]. 
Thus for a small superconducting sphere, i.e., with a ra- 
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Fig. 4: (Colour on line) Isoplots of the Cooper-pair density for 
three fi/fio = 25 accessible states, namely, 1, 2 and 3 external 
vortex pairs (left to right) in a superconducting sphere of radius 
15. 0£. Red and grey surfaces are parts of the same isosurface 
drawn at 20 % of the maximum. 



tio R/£ not very large, the presence of boundaries strongly 
influences the vortex state as seen in figs. 3 and 4. This is 
because CVLs very near to the superconductor's external 
surface eventually break and spring to the surface in the 
form of EVPs. Next we discuss some general aspects of 
these states. The switch of a CVL into an EVP induces a 
first order transition. EVPs are shown in fig. 4 in a grey 
and red color scheme. Red and grey colors describe parts 
of a single isosurface contained in distinct concentric vol- 
umes. Inside a inner cubic cell the isosurface is red and 
in the remaining volume grey. Thus a CVL is always red, 
whereas the EVPs are double colored, red and grey. The 
external isosurface, including the regions of contact with 
the EVPs are represented in grey color. Thus for the finite 
size superconductor the vortex state is generally classified 
according to the number of CVLs and of EVPs, an im- 
portant difference to previous 2D studies [12-14] where 
the only distinction is between 'giant' and 'multi' vortex 
states. States with multiple EVPs manifest as a collec- 
tion of vortices on the sample surface, but are conjoined 
to a giant vortex in the sample center. The total angu- 
lar momentum ('winding number', or 'vorticity') increases 
by one for each additional EVP. According to fig. 4 this 
giant vortex also contains three embryonic CVLs. The fol- 
lowing animations show three-dimensional views of some 
/i//io states, namely, for 24, 34, and 31. These states show 
3 embryonic CVLs with 1, 2 and 3 EVPs, respectively: 
M24.mov, M34.mov and M31.mov (see multi_media.doc 
file for more details about the animations). The CVL to 
EVP (first order) transitions lead to jumps in the magne- 
tization (M = —dF/dfi), shown in fig. 5 as a function of 
the dipole moment fi. Note that this feature directly leads 
to accessible experimental detection of EVP states, e.g. 
by Hall-probe magnetization measurements. However, the 
observation of CVL states is more demanding and requires 
a combination of various complementary microscopic ex- 
perimental techniques, such as Small Angle Neutron Scat- 
tering (SANS), muon spin rotation (/iSR) [21], and possi- 
bly scanning tunneling microscopy on curved surfaces to 
detect the threefold slight suppression of the order param- 
eter on the surface of the sphere [as shown in grey scale of 
figs. 5(a-c)]. 
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Fig. 5: (Colour on line) Sample magnetization (in arbitrary 
units) vs. dipole moment. Insets show the phase in the equa- 
torial plane of the sample for the 1-3 external vortex pair states 
of fig. 4 (for h/hq = 25, open dots), (a)-(c) 2 contourplots 
of the top half hemisphere of the sample, for the three confined 
vortex loops state and fif/Jo = 20, 24, and 28, respectively. 



The evolution of the vortex state. — The onset 
of EVPs unveils interesting aspects of the threefold sym- 
metry because of the soft breaking interactions that drive 
such first order transitions. A CVL becomes unstable near 
to the surface of the superconductor and eventually breaks 
apart giving rise to a EVP, as shown in figs. 2-5. For our 
spherical sample the three initially equivalent loops are 
naively expected to grow and reach the surface simulta- 
neously, turning into three EVPs. Nevertheless one and 
two EVPs are also observed implying that the threefold 
symmetry is disturbed during its evolution by some soft 
breaking mechanism, here caused by the use of a cubic 
grid. The evolution of the remaining two CVL, as the 
third one turns into an EVP, is quite remarkable. The 
outgrowth of one particular CVL from the threefold state 
causes: (i) the shrinkage of the other two ones and (ii) the 
appearance of a new CVL to recover the threefold sym- 
metry near the H C 2 core. This shows that it is easier to 
recover the threefold symmetry near the H C 2 core because 
there the CVLs are less energetically costly to produce. 
The development of the symmetry during the simulated 
annealing transition from to 1 EVP is shown in fig. 6. A 
3 CVL state becomes unstable at fi/fiQ = 28.5 [fig. 6(a)] 
and evolves, under controlled annealing temperature, to a 
1 EVP state. Threefold symmetry is observed both below 
and above the transition point. A trapped intermediate 
(unstable) configuration during this transition is shown 
as inset (b) of fig. 6. We confirmed these results by a 
saddle-point analysis [18], carried through the kinetic en- 
ergy eigenfunction expansion. Moreover we have shown 
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Fig. 6: (Colour on line) The free energy curves versus dipole 
moment for and 1 external vortex pair state obtained by the 
simulated annealing method. The confined vortex loop state is 
met ast able for /x//xo > 12.5 and becomes unstable at /i//io = 
28.5. The isosurfaces in the insets (a-c) show snapshots of the 
consequent to 1 external vortex pair transition. 

that the presence of an external current applied along the 
magnetic dipole can induce this transition, as shown in 
the following animation done for fi/fio = 10: ZERO-TO- 
ONE. mov (see multLmedia file for more details about the 
animations). The density increases away from the H C 2 
core at the center but eventually this growth stops and the 
trend reverts to a decreasing pattern near to the surface 
of the mesoscopic sphere. Outside the order parameter 
must rigourously vanish, and so, there is always an exter- 
nal isosurface density near to the surface of the mesoscopic 
sphere with constant and small |^| 2 . The figs. 1, 4 and 6 
show it in grey scale as part of the isosurface drawn at 20 
% of the maximum \^\ 2 . 

Conclusion. — Using several independent theoretical 
approaches, we found a threefold embryonic growth of vor- 
tex loops from the H C 2 core surrounding the magnetic mo- 
ment in a superconducting sphere. As shown here, the 
superconducting state with three vortex loops is energet- 
ically favorable over one or two vortex loops regardless of 
the sample boundaries. Therefore, our results also apply 
to bulk samples. We expect that the present scenario of 
loop breaking at the sample surface can be useful to under- 
stand properties of a bulk superconductor with some dis- 
tribution of magnetic domains in its interior, where grow- 
ing vortex loops may interconnect neighboring magnets 
rather than spring to the surface, which would lead to the 
experimentally observable spontaneous vortex phase [22]. 
The puzzling properties of the recently discovered ferro- 
magnetic superconductors [23] and superconducting ferro- 
magnets [24, 25] are indicative of vortices stemming from 
an internal magnetic field. We expect that the recent ad- 



vancements in material synthesis, which have led to several 
fabricated nano-composites, such as the MgB 2 supercon- 
ductor with embedded magnetic Fe203 nanoparticles [26], 
or Gd particles incorporated in a Nb matrix [27], can be 
used to obtain the predicted vortex patterns made of con- 
fined vortex loops and external vortex pairs. 
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Abstract. - A magnetic inclusion inside a superconductor gives rise to a fascinating complex 
of vortex loops. Our calculations, done in the framework of the Ginzburg-Landau theory, reveal 
that loops always nucleate in triplets around the magnetic core. In a mesoscopic superconducting 
sphere, the final superconducting state is characterized by those confined vortex loops and the 
ones that eventually spring to the surface of the sphere, evolving into vortex pairs piercing through 
the sample surface. 
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Examples of one- or multi-dimensional collections of 
topological defects arise in various physical systems, such 
as vortices in quantum fluids [?], dislocations in solids 
[?], global cosmic strings [?,?], and polymer chains [?]. 
Changes of symmetry in their structural and dynamical 
properties often lead to very rich phase transitions. In su- 
perconductivity, the structure and symmetry of the vortex 
lines has been a major issue even in the early days, when 
Abrikosov [?] proposed a square pattern for the lattice of 
penetrating flux tubes in bulk superconductors; later on, 
Kleiner et al. [?] found that the most efficient way to pack 
interacting parallel rods is actually through a hexagonal 
pattern. However, in systems with reduced dimension- 
ality, the Abrikosov symmetry of the vortex state can be 
strongly affected by the sample boundaries. In mesoscopic 
superconductors, where the volume to surface area ratio is 
small [?], strong lateral confinement may impose the for- 
mation of multiquanta ('giant') vortices [?,?], while the 
shape of the boundary directs the symmetry of the final 
vortex configuration [?] . 

In this Letter, we consider a novel vortex structure that 
stems from the specific inhomogeneous magnetic field of a 
point like magnetic dipole with moment fi in the center of 
a superconducting sphere. This magnet is the only source 
and sinkhole of vortex lines and sets an axis around which 
the vortex state exhibits a discrete symmetry. Because 
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the sphere does not break this degeneracy many of the 
present results also apply for a magnetic inclusion in a 
bulk superconductor. 

Therefore, the present analysis is very different from 
earlier studies of vortex and antivortex structures in su- 
perconductor/magnet hybrids [?,?,?] where magnetic dots 
were placed near, but outside a thin superconducting sam- 
ple. While in the latter case vortices were two-dimensional 
(2D) structures characterized by a coin-like core, here they 
are lines in three-dimensional (3D) space, set on a formal- 
ism that comprises the full 3D Ginzburg-Landau theory, 
thus far more complex than previous works. Consequently, 
novel vortex phases with profound 3D features are found 
for these sets of curved vortex lines, made of confined vor- 
tex loops (CVLs) and external vortex pairs (EVPs). 

For small /i, as shown later, we find that the vortex state 
is made of exactly three confined vortex loops around the 
dipole axis. For an increasing magnetic moment fi more 
elaborate vortex arrangements arise because of the inho- 
mogeneity of the field. This threefold growth of vortex 
loops in our system is shown here to be a consequence 
of the energetic balance. Threefold symmetry has been a 
source of inspiration since the sixties when its discovery 
in hadronic physics has led to quarks as the basic con- 
stituents of matter. Recently, experimental evidence was 
found for the onset of trimer bound states in ultracold 
gases of caesium atoms [?] and Z3 symmetry was also pro- 
posed as a way to solve baryon number violation in grand 
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